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^-^ In my talk I will present an overview of our recent work involving the use of supersymmetric quan- 

tum mechanics (SUSY-QM). I begin by discussing the mathematical underpinnings of SUSY-QM 
^.. and then discuss how we have used this for developing novel theoretical and numerical approaches 

d suitable for studying molecular systems. I will conclude by discussing our attempt to extend SUSY- 

^ QM to multiple dimensions. 
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I My introduction to supersymmetric quantum mechanics (SUSY) came quite 

(-H by accident. I had heard of SUSY in the context of high-energy physics where 

C^ the SUSY theory postulates that for every fermion there is boson of equal mass 

^ (i.e. energy). This comes about because for every quantum Hamiltonian there 

O^ is a partner Hamiltonian that has the same energy spectrum above the ground- 

state of the original system. In other words, above the ground state of iJi, each 
^—1 higher-lying eigenstate is partnered with an eigenstate of i?2- In particle physics, 

^ the Hi "sector" is populated by bosons and the H2 sector by fermions and SUSY 

00 predicts that the lowest lying fermion state is energetically degenerate with the 

00 first excited boson state. Evidence for SUSY has proven to be elusive and it is 

now believed that SUSY is a broken symmetry. 

Last January (Jan-09) at a conference dedicated to Bob Wyatt, my co-author 
1^ suggested that we look at SUSY as a way to develop new computational methods 

and approaches. Up until now, SUSY has been more of a mathematical technique 
that has been used more or less as a way to obtain stationary solutions to the 
Schrodinger equation for the variety of one-dimensional potential systems. In 
this paper and in my talk, I will discuss some of the work we have been doing 
in developing "SUSY" inspired methods for performing quantum many-body 



* corresponding authoi[| 



calculations and quantum scattering calculations. I shall begin with a brief 
overview of the SUSY theory and some of its elementary results. I shall then 
discuss how we have used the approach to develop both analytical and numerical 
solutions of the stationary Schrodinger equation. I will conclude by discussing 
our recent extension of SUSY to higher dimensions and for scattering theory. 



II. MATHEMATICAL CONSIDERATIONS 

Before discussing some of our recent results, it is important to introduce briefly 
the mathematical formulation of SUSY quantum mechanics. 

A. Hamiltonian formulation of SUSY 

In quantum theory, there is a fundamental connection between a bound state 
and its potential. This is simple to demonstrate by writing the Schrodinger 
equation for the stationary states as 

Vi{x)~Er,^^^^dl^n=Q[^n] (1) 

where we recognize the right-hand side as the Bohm quantum potential which 
will certainly be discussed repeatedly at this conference. One of the remarkable 
consequences of this equation is that every stationary state of a given potential 
has the same functional form for its quantum potential Q. Thus, knowing any 
bound state allows a global reconstruction of the potential, V{x) up to a constant 
energy shift. 

SUSY is obtained by factoring the Schrodinger equation into the form [IHS] 

Hi; ^ A+AiP^/^ = (2) 

using the operators 

A = ^^d^ + W and A+ = ^d^ + W. (3) 

v27n v2m 

Since we can impose Aipo = 0, we can immediately write that 

W{x) = -^d,\nij,. (4) 

W{x) is the superpotential which is related to the physical potential by a Riccati 
equation. 

V(x) = W^ix) - -^W'ix). (5) 

V2m 



The SUSY factorization of the Sehrodinger equation can always be apphed in 
one-dimension. 

From this point on we label the original Hamiltonian operator and its asso- 
ciated potential, states, and energies as Hi, Vi, ipn and £"„ . One can also 
define a partner Hamiltonian, H2 = AA^ with a corresponding potential 

V2^W^ + -^W'{x). (6) 

All of this seems rather circular and pointless until one recognizes that Vi and 
its partner potential, V2, give rise to a common set of energy eigenvalues. This 
principle result of SUSY can be seen by first considering an arbitrary stationary 
solution of Hi , 

Hi^i'^=A+A^M = El^^^^\ (7) 

This implies that (Aipn ) is an eigenstate of H2 with energy £"„ since 

H^iAyjl'^) = AA+A^/i'^ = Ei'\A^i^^). (8) 

(2) 

Likewise, the Sehrodinger equation involving the partner potential H2ipn — 

(2) (2) 
En ipn implies that 

A+AA+^i^) = HiiA-^^^:^) = Ei'HA+^l^^ (9) 

This (along with £"0 — ) allows one to conclude that the eigenenergies and 
eigenfunctions of Hi and H2 are related in the following way: i?„^i ~ En , 

^i^) = -jL=A,p^^l, and i;i% ^ ^A+^(2) (10) 

for n > 0. '51] Thus, the ground state of H2 has the same energy as the first 

excited state of Hi. If this state ipo is assumed to be node-less, then ipl ex 

A'^ipo will have a single node. We can repeat this analysis and show that H2 
is partnered with another Hamiltonian, H^ whose ground state is isoenergetic 
with the first excited state of H2 and thus isoenergetic with the second excited 
state of the original Hi . This hierarchy of partners persists until all of the bound 
states of Hi are exhausted. 



B. SUSY algebra 

We can connect the two partner Hamiltonians by constructing a matrix super- 
Hamiltonian operator 



and two matrix "super-charge" operators 

Q=(yn)-A^- (12) 

and 
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A^ 




Q^=[o )=^ "+ (13) 

where a± arc 2x2 Pauh spin matrices. Using these we can re-write the SUSY 
Hamihonian as 

The operators {H, Q, Q+} form a closed algebra (termed the Witten superalge- 
bra) with 

[H,Q] = [H,Q+]=0 (15) 

{Q,Q} = {Q+,Q+} = (16) 

{Q,Q+} = H (17) 

The first algebraic relation is responsible for the degeneracy of the spectra of 
Hi and H2 and the supercharges transform an eigenstate of one sector into an 
eigenstate of the other sector. 

As an example and perhaps a better connection to the physics implied by 
this structure, consider the case of a one-dimensional particle with an internal 
spin degree of freedom and with [x,p\ = i denoting the position and momentum 
of the particle. Conserved SUSY would imply that all non-diagonal coupling 
terms between the bosonic (coordinate) and fermionic (spin) degrees of freedom 
are exactly zero. This of course is equivalent to making the Born Oppenheimer 
approximation for a two-state system coupled to a continuous field x{t). In this 
case, SUSY is preserved so long as dttp{x{t),t) = dtip{x{t),t). SUSY is broken 
when x{t)dx'4'{x{t), t) ^{) which would lift the degeneracy between the states of 
Hi and H2. 



C. Scattering in one dimension 

The SUSY approach is not limited to bound-state problems. For a one- 
dimensional scattering system, it is straightforward to apply the SUSY theory 
to determine a relation between between the transmission and reflection coef- 
ficients of the supersymmetric partners. Asymptotically, we can assume that 
W{x) — )■ W± as a; — > ±00. In the same limit, the partner potentials become 



Vi,2 — > W^. For a plane wave incident from the left with energy E scattering 
from Vi,27 we require the following asymptotic forms: 

x— y — OO 



a;— >-+oo 



We can derive a relation between the two scattering states by using the relation 

ilj'^^\k,x) = N A~^ ip2{k' , x) . For the left-hand components (x — > — oo). 

^^kx ^ JlWg'ikx ^ ^ r/ -^ _^ ^_\ gifca; _^ /^^ _^ ^_\ g-»fcxl (20) 

where in the last line we have incorporated the h/\/2m in to the normalization 
and wrote W± — W±y/2m/h. Likewise for the transmitted coefficients (x — > 

+00). 

XW^ik'x ^ ^(_jfc' ^ ^^)y(2)g»fe'a; (21) 

Eliminating the common normalization factor and using the fact that k = 
y/2m{E — W-)/h and k' = yj2ra{E — W+)/h from the Schrodingcr equation 
we can arrive at 

TW(fc)=.^^^r(^)(fc). (23) 

W^ — ik 

Consequently, knowledge of the scattering states of Vi allows one to easily con- 
struct scattering states for the partner potential. 



D. Non-stationary states 

Finally, one can use the SUSY approach in a time-dependent context by writ- 
ing 

where ■i/'''^-* is a non-stationary state in the first sector. If Vi is independent of 
time, then the superpotential must also be independent of time and so we can 
write 

ihAdtij^^^ = ihdtiAij^^^) = AA+{A^^^^) 



In other words, we have the time-dependent Schrodinger equation for the partner 
potential 

The two non-stationary states are partnered, V'^^'' oc A~^^^-^K We also note that 
these states satisfy 

and 

for some initial state ip{0)- Using the charge operators we can show that 

As above in the scattering example, one can use the dynamics of one sector to 
determine the dynamics in the other sector. 

The partnering scheme presents a powerful prescription for developing novel 
approaches for solving a wide variety of quantum mechanical problems. This 
allows one one use analytical or numerical solutions of one problem to determine 
either approximate or exact solutions to some new problem. In the sections that 
follow, I present some of our attempt to use SUSY in a numerical context. At the 
moment our numerical results are limited to one spatial dimension. As I shall 
discuss, extending SUSY to multiple dimensions has proven to be problematic. 
However, in Sec. V we present our extension using a vector-SUSY approach we 
are developing. 



III. USING SUSY TO OBTAIN EXCITATION ENERGIES AND 

EXCITED STATES 

The SUSY hierarchy also provides a useful prescription for determining the 
excited states of Hi (which may represent the physical problem of interest.) The 
first excited state of Hi is isoenergetic with the ground state of H2- Since this 
state is node-less, one can use either Ritz variational approaches or Monte Carlo 
approaches to determine this state to very high accuracy. 

Two basic tools used in computational chemistry are the Quantum Monte 
Carlo (QMC) and the Rayleigh-Ritz variational approaches. Both approaches 
yield their best and most accurate results for ground state energies and wave 
functions. Although the variational method also gives bounds for the excited 
state energies as well as the ground state (the Hylleraas-Undheim theorem ^), 
it is well known that their accuracy is significantly lower than that of the ground 



state. Even more serious, the wave functions are known to converge much more 
slowly than the energies. 

In the case of the QMCjSHin]. there are additional difficulties associated with 
the presence of nodes in the excited state wave functions 11 . While some 
progress has been made in dealing with this issue (e.g., the "fixed node" or 
"guide wave" techniques) [8HT2] the computational effort required is greater and 
the accuracy is lower and in fact, no general solution to the diflficulty has been 
found for reducing the computational effort and increasing the accuracy for 
excited state calculations in QMC to the same level as is attained for the ground 
state. In fact, it is very likely the presence and effects of nodes in the excited 
states that is largely responsible for the lower accuracy and slower convergence 
of excited state results in the variational method. The precise determination 
of nodal surfaces is expected to play a crucial role since they reflect changes in 
the relative phase of the wave function. Because of the ubiquitous importance 
of both the variational and QMC methods, solving the so-called "node problem 
will have enormous impact on computational chemistry. 



A. Using SUSY to improve quality of variational calculations 

We now turn to the proof of principle for this approach as a computational 
scheme to obtain improved excited state energies and wave functions in the 
Rayleigh-Ritz variational method. We should note that these results can be 
generalized to any system where a hierarchy of Hamiltonians can be generated 
because of the nature of the Rayleigh-Ritz scheme. In the standard approach one 
calculates the energies and wave functions variationally, relying on the Hylleraas- 
Undheim theorem for convergence [3]. This, however, is unattractive for higher 
energy states because they require a much larger basis to converge to the same 
error. We stress that this is true regardless of the specific basis set used. Of 
course, some bases will be more efficient than others but it is generally true that 
for a given basis, the Rayleigh-Ritz result is less accurate for excited states. We 
address this situation by solving for ground states in the variational part of the 
problem. 

To demonstrate our computational scheme, we investigate the first example 
system from the previous section. For the potential 

Vi [x) :^ x^ + ix^ + x^ - 2. (24) 

exact solutions are known for all states of Hi. We use the exact results to 
assess the accuracy of the variational calculations. Here we employed a n-point 
discrete variable representation (DVR) based upon the Tchebchev polynomials 
to compute the eigenspectra of the first and second sectors. [T31 [Ti] In Fig. [T]we 
show the numerical error in the first excitation energy by comparing El (n) and 
Eq^u) from an n point DVR to the numerically "exact" value corresponding to 
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FIG. 1: Convergence of first excitation energy El for model potential Vi = x^ + 
4a;^ + x^ — 2 using a n-point discrete variable representation (DVR). Gray squares: 
e = logjo lEiin) — El{exact)\, Black squares: e = logj,, \EQ{n) — El{exact)\. Dashed 
lines are linear fits. (From Ref[15j.l 



a 100 point DVR, 



Likewise, 



elin) ^\og^o \Elin) - El{exact)\. 



el{n) = logio \El{n) - El{exact)\. 

For any given basis size, Eq < e\. Moreover, over a range of 15 < n < 40 
points, the excitation energy computed using the second sector's ground state is 
between 10 and 100 times more accurate than El{n). This effectively reiterates 
our point that by using the SUSY hierarchy, one can systematically improve 
upon the accuracy of a given variational calculation. 



B. Monte Carlo SUSY 

Having defined the basic terms of SUSY quantum mechanics, let us presume 
that one can determine an accurate approximation to the ground state density 
po (x) of Hamiltonian Hi . One can then use this to determine the superpotential 
using the Riccati transform 



W!:'^ 



h 9 In Po 
/2to dx 



(1) 



(25) 



and the partner potential 

fc2 q2i (1) 

Certainly, our ability to compute the energy of the ground state of the partner 
potential V2 depends on having first obtained an accurate estimate of the ground- 
state density associated with the original Vi. 

For this we turn to an adaptive Monte Carlo-like approach developed by Mad- 
dox and Bittner.[16j Here, we assume we can write the trial density as a sum 
over N Gaussian approximate functions 

PT(x) = ^G„(a;,c„). (27) 

n 

parameterized by their amplitude, center, and width. 

G„(x,{c„})=c„oe-^-(^-^"^)' (28) 

This trial density then is used to compute the energy 

E[pt] - (Vi) + {Q[pt]) (29) 

where Q[pt] is the Bohm quantum potential, 

2m y/px ox'^ 



The energy average is computed by sampling pt{x) over a set of trial points 
{xi\ and then moving the trial points along the conjugate gradient of 

E{x) = Vi{x)+Q[pt]{x). (31) 

After each conjugate gradient step, a new set of €„ coefficients are determined 
according to an expectation maximization criteria such that the new trial density 
provides the best A'^-Gaussian approximation to the actual probability distribu- 
tion function sampled by the new set of trial points. The procedure is repeated 
until 5{E) = 0. In doing so, we simultaneously minimize the energy and opti- 
mize the trial function. Since the ground state is assumed to be node-less, we 
will not encounter the singularities and numerical instabilities associated with 
other Bohmian equations of motion based approaches. [16H2T] Moreover, the 
approach has been extended to very high-dimensions and to finite temperature 
by Derrickson and Bittner in their studies of the structure and thermodynamics 
of rare gas clusters with up to 130 atoms. [22l |23] 
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IV. TEST CASE: TUNNELING IN A DOUBLE WELL POTENTIAL 

As a non-trivial test case, consider the tunneling of a particle between two 
minima of a symmetric double potential well. One can estimate the tunneling 
splitting using semi-classical techniques by assuming that the ground and excited 
states are given by the approximate form 

i^± = ^iM^)±M-x)) (32) 

where (JJq is the lowest energy state in the right-hand well in the limit the wells 
are infinitely far apart. From this, one can easily estimate the splitting as |24j 

6^4-UOmO) (33) 

m 

If we assume the localized states {(t>o) to be gaussian, then 

^± ex ^(e-/5(— ")' ± e-/5(-+-°)') (34) 

and we can write the superpotential as 

W = \ —h^ (x - Xo tanh(22;Xo;3)) . (35) 

V m 

From this, one can easily determine both the original potential and the partner 
potential as 

^1,2 = W^± -^W (36) 

v2to 



m 

„2 



2{x — Xo tanh(2a;Xo/3)) 



2 



± (2x^sech^(2a;Xo/3) - 1) (37) 

While the Vi potential has the characteristic double minima giving rise to a 
tunneling doublet, the SUSY partner potential V2 has a central dimple which 
in the limit of Xq —>■ cxs becomes a J-function which produces an unpaired and 
node- less ground state. |3] Using Eq. M one obtains ipi — ip^ (x A^^o which 
now has a single node at a; = 0. 

For a computational example, we take the double well potential to be of the 
form 

Vi{x)=ax^ + bx^ + Eo. (38) 
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FIG. 2: (a) Model double well potential(blue) and partner potential (purple). The 
energies of the tunneling doublets are indicated by the horizontal lines at V — Ocm"'^ 
and V = 59.32 cm~^ indicating the positions of the sub-barrier tunneling doublet, 
(b) Final ground state density (blue) superimposed over the Gaussians used in its 
expansion, (purple) (From Ref. ,25j) 



with a = 438.9cm-^ /{bohr^), b = 877. 8cm-^ /{bohr)'^, and Eo = -181.1cm-i 
which (for m = mn ) gives rise to exactly two states at below the barrier 
separating the two minima with a tunneling splitting of 59.32 cm~^ as computed 
using a discrete variable representation (DVR) approach. [55] For the calculations 
reported here, we used Up = 1000 sample points and iV = 15 Gaussians and in 
the expansion of pt{x) to converge the ground state. This converged the ground 
state to 1 : 10^* in terms of the energy. This is certainly a bit of an overkill in 
the number of points and number of gaussians since far fewer DVR points were 
required to achieve comparable accuracy (and a manifold of excited states) . The 
nunrerical results, however, are encouraging since the accuracy of generic Monte 
Carlo evaluation would be Xj Jrhp ~ 3% in terms of the energy. [52^ Plots of V\ 
and the converged ground state is shown in [2] 

The partner potential V2 = W"^ + hW' j-^phn^ can be constructed once we 
know the superpotential, W{x). Here, we require an accurate evaluation of the 
ground state density and its first two log-derivatives. The advantage of our 
computational scheme is that one can evaluate these analytically for a given set 
of coefficients. In[2k we show the partner potential derived from the ground-state 
density. Where as the original V\ potential exhibits the double well structure 
with minima near Zq = ±1 , the V2 partner potential has a pronounced dip about 
a: = 0. Consequently, its ground-state should have a simple "gaussian"-likc form 
peaked about the origin. 

Once we determined an accurate representation of the partner potential, it is 
now a trivial matter to re-introduce the partner potential into the optimization 
routines. The ground state converges easily and is shown in [3k along with 
its gaussians. After 1000 CG steps, the converged energy is within 0.1% of 
the exact tunneling splitting for this model system. Again, this is an order of 
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FIG. 3: (a) Ground state density of the partner Hamiltonian H2 (blue) superimposed 



(1) 



over its individual Gaussian components, (b) Excited state ip\ derived from the 
ground state of the partner potential, tpo ■ (From Ref . 25] ) 

node position 




CO Step 
100 200 300 400 500 600 ^ 

FIG. 4: Location of excited state node for the last 600 CG steps. (From Ref. [25)) 



magnitude better than the Xj Jn^ error associated with a simple Monte Carlo 
sampHng. Furthermore, Wa shows tA^ 



A^^-00 computed using the converged 



\ density. As anticipated, it shows the proper symmetry and nodal position. 
By symmetry, one expects the node to lie precisely at the origin. However, 
since we have not imposed any symmetry restriction or bias on our numerical 
method, the position of the node provides a sensitive test of the convergence 
of the trial density for pp . In the example shown in FigM the location of 
the node oscillates about the origin and appears to converge exponentially with 
number of CG steps. This is remarkably good considering that this is ultimately 
determined by the quality of the 3rd and 4th derivatives of ph that appear 
when computing the conjugate gradient of V2. We have tested this approach on 
a number of other one-dimensional bound-state problems with similar success. 
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V. EXTENSION OF SUSY TO MULTIPLE DIMENSIONS. 

While SUSY-QM has also been explored for one dimensional, non-relativistic 
quantum mechanical problems [51 [?7H5T] . thus far these studies have focused on 
the formal aspects and on obtaining exact, analytical solutions for the ground 
state for specific classes of problems. In several recent papers[T3[151[311[33], we 
have begun exploring the SUSY-QM approach as the basis of a general computa- 
tional scheme for bound state problems. Our initial studies have been restricted 
to one dimensional systems (for which there are, obviously, many powerful com- 
putational methods). In our first paper, we found that SUSY-QM (combined 
with a new periodic version of the Heisenberg-Weyl algebra) yields a robust, 
natural way to treat an infinite family of hindered rotors. [55] Next we showed 
that the SUSY-QM leads to a general treatment of an infinite family of an- 
harmonic oscillators, such that highly accurate excited state energies and wave 
functions could be obtained variationally using significantly smaller basis sets 
than a traditional variational approach requires[15j. Most recently , we have 
considered a 1-D double well potential in which we solved for the ground state 
energy and wave function using a VQMC approach. Then using SUSY-QM, we 
(numerically) generated an auxiliary Hamiltonian whose nodeless ground state 
is iso-spectral (degenerate) with the first excited state of the original system 
Hamiltonian. This ground state was also easily determined by VQMC, yielding 
excellent accuracy for the first excited state energy. |25j Even more significant, 
by using the charge operators naturally generated in the SUSY-QM approach, 
we also obtained excellent accuracy for the first excited state wave function. 
Furthermore, at no point did impose a fixed node or symmetry on the excited 
state wave function and our calculation only involved working with a nodeless 
ground state. 

Of course, all this begs the question: Can this approach be generalized to 
higher numbers of dimensions and to more than a single particle? There has 
been substantial effort in the past to do just this.^, 'SS'-'W, '34'-^45; However, 
to date, no such generalization has been found that is able to generate all the 
excited states and energies even for so simple a system as a pair of separable, 
1-D harmonic oscillators (HO) or equivalently, for a separable 2-D single HO. 
In our most recent, unpublished work [32], we have succeeded in obtaining such 
a generalization and showed that it does, in fact, yield the correct analytical 
results for separable and non-separable problems. In the next section, we present 
a succinct summary of our approach. The major question now is whether this 
formalism provides a basis for a robust, computational method for determining 
excited state energies and wave functions for large, strongly correlated systems 
using either QMC or variational algorithms applied solely to nodeless ground 
state problems. 
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A. Difficulties in extending beyond one dimension 

To move beyond one dimensional SUSY, loffe and coworkers have explored 
the use of higher-order charge operators [HI 021 HHl HZ] , and Kravchenko has 
explored the use of Clifford algebras 48 . Unfortunately, this is difficult to do in 
general. The reason being that the Riccati factorization of the onc-dimcnsional 
Schrodinger equation does not extend easily to higher dimensions. One remedy 
is write the charge operators as vectors A = {+d + W) and with A^ = {—d + 
Wy as the adjoint charge operator. The original Schrodinger operator is then 
constructed as an inner-product 

Hi= A+ ■ A. (39) 

Working through the vector product produces the Schrodinger equation 

Hic/} = (-V^ + W^ -{V ■ W))(j) = (40) 

and a Riccati equation of the form 

U{x) = W^ ~V -W. (41) 

For a 2d harmonic oscillator, we would obtain a vector superpotential of the 
form 

W^ \-Wi^l^'>^{x,y)^{W^,Wy) (42) 

Let us look more closely at the V • W part. If we use the form that W = — V In i/i, 
then —V • Vln?/; = — V^ \nijj which for the 2D oscillator results in V • W" = 2. 
Thus, 

W'^ -^ -W ^{x^ +y^)-2 (43) 

which agrees with the original symmetric harmonic potential. Now, we write 
the scaled partner potential as 

U2 = W^ +V -W ^{x^ + y^) + 2. (44) 

This is equivalent to the original potential shifted by a constant amount. 

U2 = Ui+ 4. (45) 

The ground state in this potential would be have the same energy as the states 
of the original potential with quantum numbers n -{- m = 2. Consequently, even 
with the this naive factorization, one can in principle obtain excitation energies 



15 

for higher dimensional systems, but there is no assurance that one can reproduce 
the entire spectrum of states. 

The problem lies in the fact that neither Hamiltonian H2 nor its associated 
potential U2 is given correctly by the form implied by Eq. |40| and Eq. |44[ Rather, 
the correct approach is to write the H2 Hamiltonian as a tensor by taking the 
outer product of the charges H2 = AA^ rather than as a scalar A ■ A^. At 
first this seems unwieldy and unlikely to lead anywhere since the wave function 
solutions of 

H2if = E^ (46) 

are now vectors rather than scalers. However, rather than adding an undue com- 
plexity to the problem, it actually simplifies matters considerably. As we demon- 
strate in a forthcoming paper, this tensor factorization preserves the SUSY alge- 
braic structure and produces excitation energies for any n— dimensional SUSY 
system. Moreover, this produces a scalar 1— >■ tensor i-> scalar hierarchy as one 
moves to higher excitations. [32] 



B. Vector SUSY 

We now give a brief summary of our new generalization of SUSY-QM to treat 
higher dimensionality and more than one particle. Previous attempts generally 
involved introducing additional, "spin-like" degrees of freedom . [^ [5U1 157^01 1^ 
H51 HH HZl SH] In our approach, we make use of a vectorial technique that can 
deal simultaneously with either higher dimensions or more than one particle. In 
fact, the two problems are dealt with in exactly the same manner. Therefore, for 
simplicity, we consider a general n-dimensional distinguishable particle system 
with orthogonal coordinates {a;^}. The Hamiltonian is given by |53| 

H = -V^ + Voix,,---Xn) (47) 

and the nodeless ground state satisfies the Schrodinger equation, 

h4'^=E^,'^4'\ (48) 

We now define a "vector super-potential" , Wi , with components 

Wi^^-—ln4'\ (49) 

Then it is easily seen that the original Hamiltonian can be recast as 

Hi^{-^ + Wi)-{^ + Wi) = Q+ -Qi (50) 
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where the Qi and Q{ are niuhi-diniensional generahzations of the SUSY charge 
operators from Eq. [2] This defines our "sector-1" (or "boson") Hamihonian 



and Eq. ( 48 1 can be written as 



FiV-r - 4'V^'^ (51) 

One can show that the vector superpotential is related to the original (scalar) 
potential via: 

Vo = Wi-Wi-V-Wi. (52) 

The various components of the charge operators, Ai and Af are defined by 

A,, ^A_^+W,,kAl = -£-^+ W,,. (53) 

Note that since these are associated with orthogonal degrees of freedom, the 
charge operators can be applied either by individual components or in vector 
form. 

Next, consider the Schrodinger equation for the first excited state of H. We 
can write this using the charge operators as 

H,4'^ = Ei'^4'^ = {A+ ■ A, + 4^ Vr^ (54) 

We apply Ai to Equation (9): 

iA,A+) ■ A.^p^"^ = (i?W _ i?("))IiV('^ (55) 

Here we identify (Ai^j) as a new, auxiliary Hamiltonian. It is important to 
note that this is constructed from the outer or tensor product of the charge 
operators rather than from inner or dot product as used in constructing Hi. 

Its eigenvector, Ai^pi , is isospectral with the excited state, -0^ of i7i ( since 

Eq is known, determining {E[ — E[ ) yields E^ ). We therefore define the 
tensor Hamiltonian for the second sector as 

1/2 = AiA\ (56) 

and vector state function as 

It is easy to show that the ground state energy of ii 2 is related to the first 
excitation energy of the original Hamiltonianm 

E^^'^ ^ e['^ - E^fl (58) 
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Furthermore, the ground state of H 2 is also nodeless. This has been exphcitly 
shown to be true for the separable 2-particle HOs considered earlier[32]- There- 
fore, we propose to apply both the VQMC and the standard variational methods 
to determine E^ and V'o • Of course, knowing the second sector ground state 
energy also gives us the first excited state energy of the original Hamiltonian 
(Eq. ISS])). Furthermore, we form the scalar product of 

^2-4'^=4'VT (59) 



with A^ obtaining 



{A+ ■ 11)1+4'^ = E'^Ut ■ 4') (60) 



Clearly, this is exactly 

HMt-^'i'') = E\'\At-r^'i^) (61) 

so we can conclude that 

41) = ^^=(1+. 4^)) (62) 

Thus we also obtain the excited state wave function without any significant 
additional computational effort. This is because applying the charge operator is 
much simpler than solving an eigenvalue problem (it is a strictly linear opera- 
tion). Evidence from our 1-D studies indicates that the accuracy of the excited 
states obtained using the SUSY-QM charge operator is significantly higher, for 
a given basis set, than what is obtained variationally (or with QMC) from the 
original Hamiltonian |151 [52] , This procedure can be continued as follows. We 
define a sector-2 vector super-potential with components 

(63) 



(64) 





^^^-Al^^<' 




Then it follows that 








M ■ 4'^ = (V + W2) ■ 4'^ 


= 


so we can write 


^2^A+A2+E'ih 





(65) 



and Eq. ( 59 ) is still satisfied. We form the scalar product of A2 with the first 



excited state Schrodinger equation to obtain 

(I2 • A+)A2 ■ 4"' = EfU2 ■ 4'^ (66) 
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Then we define the sector 3 scalar Hamihonian by 



H3=A2-A++ 4') (67) 



with the ground state wave equation 

i?34'^ - 4'V^'^ (68) 

It is easily seen that Eq = El — Eq. This procedure continues until all bound 
states of the original Hamiltonian are exhausted. It should also be clear that 
the sector 2 excited state wave function is obtained from the nodeless sector 3 
ground state by applying ^J to it. Then the second excited state for sector 
1 results from taking the scalar product of Af with ipi ' . The approach thus 
leads to an alternating sequence of scalar and tensor Hamiltonians, but in all 
cases we need only determine nodeless ground states. 

There are two additional aspects of the tensor sector problem that require dis- 
cussion. First we consider the validity of the Rayleigh-Ritz variational principle. 
It is easily seen from Eq. (56) that H 2 is a, Hermitian operator. Therefore, 



its eigenspectrum is real and its eigenvectors are complete. With these facts in 
hand, the proof of the variational principle follows the standard one in every 
detail. This is also true for the Hylleraas-Undheim theorem. 

Second, the QMC method is also directly applicable to the tensor sector prob- 
lem. For the example discussed above, we note that the energy is given by 

/ driftrial ■ ff 2 • '4'trial ,„„x 

Etrial = TT-r-^ f (oy) 

J drtptrial ■ Vtrial 

We next note that the integral can be expanded in terms of its components as 

It is then clear that each separate integral can be evaluated by QMC. For ex- 
ample, the [i^ V cross term is divided and multiplied by 



ll^p,,trial / dT'll)^^trial'^fi..trial- (71) 

Then the sampling is done relative to the mixed probability distribution, 

T-j W 11. trialWu. trial /'70^ 



A similar expression applies to each term in the energy expression and the eval- 
uation would need to be performed self-consistantly. 
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Thus far, we have developed a formahsm that appears to be suitable for ex- 
tending the SUSY-QM technique to higher dimensional systems. We believe 
the approach we have outlined above will provide the mathematical basis for 
a number of potentially interesting theoretical results. Moreover, we anticipate 
that when combined with either variational or Monte Carlo methods, our multi- 
dimensional extension of SUSY-QM will facilitate the calculation of accurate 
excitation energies and excited state wave functions. 



VI. OUTLOOK 

I presented a number of avenues we are actively pursuing with the goal of 
using SUSY-QM or SUSY-inspired-QM to solve problems that are difficult to 
solve using more conventional approaches. In addition to what I have discussed 
here we exploring the use of the Riccati equation to solve quantum scattering 
problems. It is as if one of the co-authors of this paper (DJK) has come full- 
circle since one of his first papers concerned solving the Hamilton-Jacobi for the 
action integral in quantum scattering. |50j [55] 



iS/h = - / W{r')dr'. 



The integrand in this last equation is the SUSY super-potential. Furthermore, 
there is a connection between our work and the complex-valued quantum tra- 
jectories studied by Wyatt and Tannor and their respective co-workers. 
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